Effect of Induced Spin-Orbit Coupling for Atoms via Laser Fields 
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We propose an experimental scheme to observe spin-orbit coupling effects of a two-dimensional 
(2D) Fermi atomic gas cloud by coupling its internal electronic states (pseudospins) to radiation 
in a Lambda configuration. The induced spin-orbit (SO) coupling can be of the Dresselhaus and 
Rashba type with and without a Zeeman term. We show that the optically induced SO coupling 
can lead to a spin-dependent effective mass under appropriate condition, with one of them able to 
be tuned between positive and negative effective masses. As a direct observable we show that in 
the expansion dynamics of the atomic cloud the initial atomic cloud splits into two clouds for the 
positive effective mass case regime, and into four clouds for the negative effective mass regime. 

PACS numbers: 71.70.Ej, 37.10.Vz, 03.75.Ss, 05.30.Fk 



Spin-orbit (SO) coupling effect in semiconductors has 
emerged in the solid-state community as a very active 
field of research, fuelled in part by the field of spin- 
tronics [H, e.g. the engineering of devices where the 
spin degree of freedom of the electron is exploited for 
improved functionality. This has lead to new develop- 
ments in the anomalous Hall effect (AHE) 0] and the 
spin Hall effect (SHE) 0, i]. In correspondence to the 
spin of an electron, the internal degree of freedom of an 
atom (pseudospin) is represented by the superposition 
of its electronic states (hyperfine levels). SO couphng 
can be equivalently depicted as the interaction between 
an effective non-Abelian gauge potential and a particle 
with (pseudo)spin. In quantum systems, the idea gener- 
ating a gauge field adiabatically was proposed by Wilczek 
and Zee more than twenty years ago |5|. Recently, such 
an idea was applied to atomic systems, where the mo- 
tion of atoms in a position dependent laser configura- 
tion gives rise to an effective non-Abelian gauge poten- 
tial [1, 0, H, H which can lead to an effective SO 
interaction in an ultracold atomic gas 11 , 13, l3| • 



Realization of SO interaction in atomic gases opens 
new possibility of studying spintronic effects, e.g. spin 
relaxation [lH, Zitterbewegung [l3| and SHE, in atomic 
systems which provide an extremely clean environment, 
allowing in a controllable fashion unique access to the 
study of complex physics. However, experimental detec- 
tion of such SO effects in atoms requires to measure the 
pseudospins (not just hyperfine levels) that are usually 
not directly observable for atomic systems. In this let- 
ter, we propose an experimental scheme to study SO cou- 
pling effects, based on a trapped two-dimensional (2D) 
Fermi atomic gas with a simple internal three-level A- 
type setup. We demonstrate that an effective SO inter- 
action, e.g. Rashba and linear Dresselhaus terms, can 
be obtained by coupling atoms with a three-level con- 
figuration to spatially varying laser fields. The optically 
induced SO coupling can lead to a spin-dependent effec- 
tive masses under proper condition. A direct observable 
of this effects is in the expansion dynamics for each of the 



effective mass cases after the external trap is switched off 
and we predict that the initial atomic cloud splits into 
two or four clouds. 
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FIG. 1: (Color online) (a) Three-level A-type system cou- 
pled to position-dependent laser fields with large detuning; 
(b) Laser configuration for fii. 

We consider a cloud of quasi 2D {y-z plane) Fermi 
atomic gas with internal three-level A-type configuration 
(see Fig. 1 (a)) coupled to radiation. The transition 
1 6) — > I a) is coupled by the laser field with Rabi-frequency 
rii = Slio exp[i(/<i(r)] and the transition |c) —>■ \a) is cou- 
pled by another laser field = O20 exp[i02(r)], where 
(/)i_2(r) are position-dependent phases. The Hamiltonian 
of a single particle reads: H = Hq + I4rap(r) + Hj, where 
Vtrap(r) is the 2D harmonic trap, and the interacting 
Hamiltonian is given by 

Hi = hA\a){a\ - {hni\a){b\ + hn2\a){c\ + h.c). (1) 

Diagonalizing this Hamiltonian yields three eigenstates: 
\xd) = sm6\b) — cos0e"^|c), \xbi) = cos a cos 6*16) -I- 
cosasin0e*'^|c)-|-sinQ;e"^i|a), and IXB2) = sinacos0|6)-|- 
sin a sin 0e"^|c) — cosae*'^^ |a). Here cj) = 4>i — (j>2, 
the mixing angles tan^ = r22o/^io f^nd tana = 
\/^Iq + f^2o/^ = ^o/A- The corresponding eigenvalues 
are Ed = and Eb,,^ = ft(A T y/A^+Anl)/2. Since 
spatially-varying lasers arc employed, diagonalization of 
the interactirig Hamiltonian Hj leads to a SU(3) gauge 
potential 0, 0, S, We consider the large detuning 
case, > fig, where \Ed - Eb^ \ < ^0 and \xd), IxBi) 



2 



spans a near-degenerate subspace, with their eigenval- 
ues far separated from that of Eb2- We can then ap- 
ply the adiabatic condition by neglecting the state \xb2)^ 
which leads to a U(2) non-Abelian adiabatic gauge po- 
tential based on the near-degenerate subspace spanned 
by IX-D.-Bi) [14 1 . T his situation is different from the cases 
in Refs. |9| and I 111 where the adiabatic condition is also as- 
sumed between the states \xd) and Ixsi): and thus the 
spin-dependent gauge potential is still Abelian in their 
case The adiabatic non-Ablelian gauge potential for 
the present near-degenerate subspace is obtained via 

he 

A(r) - i — {xD\ ® (xbiIVIxbJ ® |XD>. (2) 
e 

For our purpose, we shall set the parameters (f>i — (1)2 ~ 
kiy and 6 — k2Z, which means fli^2 are standing waves in 
the z direction but plane waves in the y direction. Such 
configuration can be achieved by applying two laser fields 
for each atomic transition. For fli, for instance, we can 
set two laser fields with the same strength, where one 
travels with the wave vector kiCy + fc2ez and the other 
travels with kiCy — k2ez (see Fig. 1 (b)), similarly can 
be done for The Rabi- frequencies are then given by 
ni{r) = nocosik2z)e'''^y and f72(r) = nosm{k2z)e'^^y . 
For simplicity, in what follows, we use the spin lan- 
guage and denote by | t) = \xd), \ i) = Ixsi)- Under 
this configuration the gauge field ([2]) can be recast into 
A = mXiUzey — m\20yez ~ mXiIiy with the coefficients 
Ai = ?ifcifio/(2mAo) and A2 = hk2/m (note A does not 
depend on mass m, and its current form is for the defini- 
tion of Ai_2)- In addition, the scalar potentials are given 



£^^kl and ^xT = 0. 



To this step we can obtain the effective Hamiltonian 
for the near-degenerate subspace with SO coupling in the 
form (neglecting constant terms): 



H = Hn 



Hz + Vt 



rap 1 



(3) 



p2 p2 

where = + 1^, H^o = -Mc^zPy + >^2<yyPz, 

MqUz with Mq = ^-^k2 -t- ft2A' 2D harmonic 

trap Vtrap = ^mLj^{y^ + z^). A Hamiltonian of this form 
is predicted to give SHE [3| for Mq = and planar Hall 
effect [15[ for Mq ^ in solid state systems. 

The term Hgo + Hz can be readily diagonal- 
ized in the momentum space, {Hgo + HzY^°'^ = 



XikyY 



Xlkloz 



where fcj, ^ are wave vectors 



of atoms in the y and z directions. The associated eigen- 
states are given by |-|-) — [cos i?/2, i sin?9/2]^, |— ) = 
sin i?/2, cos 1^/2]^ with tani? = X2hkz/{Mo - Xihky). 
However, even in this situation, the full Hamiltonian H 
is not diagonalizable. Practically, we can consider the 
case where > Xlh^ik'pf > Xlh^ik^pf (the vaHd- 
ity in realizable experimental set-ups will be discussed 
below), kp^ are the y/z-components of the Fermi mo- 
menta, and expand the term (TJ^o + HzY^"^^ to the fc^ 



order: {Hso + Hzf'^^ ~ (Mq 
tuting this result into eq. ([3]) yields 
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(4) 



m ;^ + ^cr^. Eq. dH) shows that the SO cou- 



where -i 

777 

pling leads to a spin-dependent effective mass of atoms 
moving in the z direction, with m± = (1 ± Sm/'m)~^m 
and Sm — ?i^fc|/Mo. The eigenf unction of Eq. (j4|) reads 

K,,,Jy,z) = ^nAujy,y)rnA^2,^)e''^^'"''^^''\v), (5) 

with I77) = |±). Here 'tpuyjipn^ harmonic oscillator 
wave functions. The eigenvalues are given by E^l^ „^ — 
{uy + ^)hhjy + {uz + ^)huj^ + rjAIo — ^mXl, where — uj, 
= ijjyjl + rjSm/m are spin-dependent effective trap 
frequencies due to SO coupling, and ny,nz are integers. 

The results (jlj and (O are valid for the situation 
dm < m, where the effective mass of particles in state 
|— ) is positive, i.e. m_ > 0. On the other hand, the 
negative effective mass regime for the |— ) state can be 
reached when Sm > m. In this case one can verify the 
dispersion relation e'/^ = Hq + Hso'^ + Hi~^ (for atoms 
in |— )) represents a double- well potential in the kz axis, 
therefore a higher-order expansion with momentum kz is 
required to derive the effective Hamiltonian. Similarly, 
we consider that X^K^ikpY ^ and can then expand 
EjT up to the k^ order, which gives a double well 0'*- 
type potential form in fc^-momentum space (to_ < 0), 



I.e. 



8M3 



-I- 2^Pz 'l^l. Higher-order terms in 
the expansion can equivalently lead to small corrections 
to the coefficients of fc^ and fcf terms. Finally we get 
effective Hamiltonian for atoms at state I—): 



H 



(2) 



1 

2m 
1 



2m- 



{Py + Ximf 
-(l+72)n' 



2L 
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{l+li)Pt 



-mu'^r'^ - imA?(6) 
2 2 ' 



Here the coefficient correction 71 — 



1,72 



^ I '^Sm — 1 are small under the condition A|fi,'*(fc|-)'* <C Mq, 

say, 7^ 2 1 in the present case. 

Unlike the positive mass regime, calculation of the 
Fermi energy in the negative mass case cannot be done 
analytically. Nevertheless, we are interested in systems 
with a large number of atoms, where the Thomas-Fermi 
(TF) approximation is suitable [13, 18, [l^l- Then the 
Fermi energy can be calculated by solving the equation: 



N = 



k<kl,r<R'l 



dWkp,f (r,k,r = 0), (7) 



where the atomic distribution function in phase space 
is given by p^'' = (27r)-2(e'^(^±'('"''')-^^) + 1)"! with 
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f3 = l/ksT, and the initial size of the atomic cloud is 
given by i?| = (2ep. T 2Mo - 2miii{e'^ })^^'^ /{m^/^uj). 
From Eq. ([7]) the relation of Fermi energy to the number 
of atoms N, trap frequency and SO coupling strength 
can be obtained numerically as shown in FigO 
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|k,c) « i|k + /c2e2,-) + i|k-fc2e^,-)-^|k + A:2e,,+ 
^|k — ^26^, +). We find that the atomic density is: 

-}2 



rib 



(r,r,t) = E / ^|p^^'^[r--(k + 7?M.),k,r,o] 

+ p(^')[r--(k-r;fc2e.),k,T,0]l, (8) 
' m J 



where p^' (i"; k, T, 0) denote the distribution functions for 
atoms in the state |±) for the positive mass (j = 1) and 
negative mass (j = 2) cases. The function nc{r,T,t) 
can be calculated in the same way. Practically, we can 
assume that before the expansion begins ep < Mq, and 
then p^j^^ (r, k, r, 0) = 0. For the positive mass case, the 
evolution of the atomic density can be calculated exactly: 



FIG. 2: (Color online) Number of atoms versus Fermi en- 
ergy at 5m/m — 1.07, 1.13, 1.20, corresponding to the effec- 
tive mass m_ — —15m, —7.7m, —5m. 

There are specific observables of this spin-dependent 
effective mass induced by the spin-orbit coupling. Firstly, 
the anisotropy of the effective mass can lead to the 
anisotropic momentum distribution. Considering the 
Thomas-Fermi approximation, at zero temperature and 
for the positive mass case, the momentum distribution is 
given by n±(k) = ^[sj.TMo-£,{kl+£^kl)]. This 
contrasts with the result for the usual isotropic mass, 
where the atomic distribution in the momentum space 
is always isotropic irrespective of the shape of the trap 
17|. For the negative effective mass case, n+(k) has the 
same form, while n_(k) — „ ^ i lep + Mn — ^ik"^, + 

;f^(H-72)fc2-f ^^(l-h7i)fc^)]. Also, a fully polarized 
Fermi gas is obtained when ep < Mq, and n+(k) = 0. 
The anisotropy of the momentum distribution can lead 
to anisotropy of the Fermi velocity in the y-z plane. For 
instance, when 5m /m = 3/4 and ep < Mq, we have 
Vp « The anisotropy of Fermi velocity can be 

directly detected by time-of- flight absorption 18. loj. 

However, the dramatic signature of SO effects in the 
present Fermi atomic gas resides in the expansion dynam- 
ics of the atomic cloud after the 2D external trap and the 
laser fields are switched off. The evolution of the atomic 
distribution in phase space p(r, k, T, t) can be calculated 
by the Boltzmann transport equation assuming that at 
t = 0, Vtrap and Mq 0. For the present cold dilute 
non-interacting Fermi gas, the evolution of Pb,c{'r, k, T, t) 
is followed by the ballistic law /9{,_c(r, k, T, i > 0) = 
Ph,c(r — Ski/m, k, T, 0) (here |k, = 6, c) are the 

eigenstates after the trap and laser fields are turned 
off). The temporal atomic spatial density can be cal- 
culated by: nb,c(r,r, t) = / d'^\s.pb^c(j — Kkt/m,'k,T,0). 
The eigenstates |k, a), {a — b,c) are related to the ini- 
tial pseudospin basis |±) by |k, 6) w i|k -I- k^ez-.V) + 



Ik - k2ez,+) 



fee 



5j|k - fee^, -), and 



nb{r,T,t) 



^JrwmZ 



2j-2 



,1/2 



(9) 
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-/3iJ+(r,t) 



g-/3£;+(r,t) 
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where S± (r, t) =/ij^ + Afo-2(TT^ 2(i+^..2t2) 
One can verify that at t = the maximum point of Ub 
is obtained at r = 0, while after a sufficiently long time 
there are two maximum points at y = 0, z = ±hk2t/m. 
As a result, Eq. Q represents an initial atomic cloud 
that splits into two clouds each moving in opposite direc- 
tion with group velocities Vg = zt^ifeez/m. Using typical 
parameters: Mq = 10^^ • s-\m w 0.963 x IQ-'^^kg (^Li 
atoms) and fe = 0.87 x lO^m"^, one finds m_ w 4m, 
and Vg « ±9.0cm/s. 

The evolution of the atomic density in the negative 
effective mass regime is more complicated. Again, as- 
suming e^- < Mq, we find from Eq. ([5]) the temporal 
atomic density as (denoting by z,, = z ± hk2t/m): 



1 + 



rib 



2m ^''y 



2mh? 



{hkz 



muj y 



1 + 



9 9 ^9 



2mm- 



2 + 2cj2t2 2m + 2m-uoH'^ 



where F = mA|fi-^(l + 7i)/(4Mg). To make a qualitative 
analysis, we first calculate the time independent momen- 
tum distribution: nb{k,T,t) = gTrmo^^ ^'^'^ ^ 



In 



l+e 



)■ Here £±(k) = fip + Mo 



-/3£+(k) 

2m 2m l^V ' 



r((fe±fc2)2-T) ] with T = -m(l-f72)/(2m_r). Note 
that if fe 7^ nb(k,T,t) has four distinct maximums 
at: kz = -fc2 - VT, -fe + VT, fe - \/T, fe + VT, 
which indicates that the initial atomic cloud is composed 
of four overlapping clouds that will travel in the z-axis 



4 



at different speeds. Using parameters similar to the pre- 
vious system: m « 0.963 x 10"^^% (^Li atoms), ^2 = 
1.0x10^1, and Mo = 0.90x lO^ft-s-i, we find that T « 
4.43 X 10-i5m2 and T = 8.9 x lO^^m-'^. Thus the four 
maximums correspond to Vg — ±13.4cm/s, ±7.4cm/s. 

Figure [3] displays numerical estimates of the splitting. 
The time at which such splitting can be observed is found 
by estimating the initial size of the atomic cloud, Rp- For 
system with TV w 10^'"^ atoms and a trap frequency lu ^ 
35HZ, Rp « 230fim in the positive effective mass case 
and Rp ~ 127.8/ito when the effective mass is negative, 
assuming T < Tp. As a result, for the positive effective 
mass regime, one can verify that the two atomic clouds 
will be fully separated after the system evolves by At ~ 
5ms. For the negative mass regime, the time is about 
At ~ 9ms. It should be emphasized that measurement 
of the present expansion dynamics needs to detect the 
density of atoms in the hyperfine level |6) (or |c)), rather 
than to detect the pseudospin states |±) that are not 
differentiable for the atomic system, thus the SO effect 
obtained here is directly observable in experiments by 
direct imaging of the separated atomic clouds. 
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FIG. 3: (Color online) Splitting of atomic cloud for the pos- 
itive effective mass case (a-c) and negative mass case (d-f) 
after turning off the trap and laser fields. The x n label rep- 
resents the value each panel has been multiplied by to keep 
figure to scale. 

Now we confirm the validity of the approximation 
M§ > Xlh^ik'p)^ > Xlh'^iklf. The second inequal- 
ity is always valid, since A2 ^ Ai in the large detun- 
ing case. Note that kp can be calculated by kp = 
y/T + 2m{ep + h^TT^ /2m)/{h'^T) for m_ < and k'p = 
y^2m-ep /h? for m,_ > . For the system composed 
of TV ~ 10^ atoms and a trap frequency u! — 35HZ, 
and using the previously employed parameters, we find 
kf, « 1.41 X lO^m-i with Xlh^ikl^f/M^ « 0.019 < 1 for 
m_ > 0, and fc|, « 3.55 x lO^m"! with A| ft" 4/^/4 ^ 
0.028 ^ 1 for m_ < (note for the case m_ < the 
approximation is up to the fc^ order). Thus, the first 
inequality is also valid. Finally we estimate the TF ap- 



proximation that is considered in the calculation. TF ap- 
proximation fails in a small peri phery region Rp — SR < 
r < Rp of the atomic cloud [l7|. For the 2D fermi atom 
gas, one can find the ratio of 6R to the atomic cloud size 
Rp satisfies SR/Rp ~ TV~^/^. Therefore such a small 
region can be safely neglected for the case with a large 
number of atoms, say N > 10*. Note this model can 
be readily extended to Bose-Einstein condensate systems 
2l|, where, together with the atom-atom interaction, the 
SO coupling may lead to intriguing new physics. 

In conclusion, we have proposed an experimental 
scheme to study SO coupling effects for a cloud of a 
2D trapped Fermi gas. Under certain conditions, the 
optically induced SO coupling in atoms leads to a spin- 
dependent effective mass which can be positive or nega- 
tive. In the expansion dynamics of the atomic cloud after 
switching off the trap, it is shown that the initial atomic 
cloud splits into two or four clouds moving oppositely de- 
pending on tunable spin-orbit coupling parameters. The 
present scheme provides an applicable way to directly 
observe the SO coupling in cold Fermi atoms. 
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